Modulation of charge-density waves by superlattice structures 
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We discuss the interplay between electronic correlations and an underlying superlattice structure 
in determining the period of charge density waves (CDW's), by considering a one-dimensional Hub- 
bard model with a repeated (non-random) pattern of repulsive (L'^ > 0) and free (L^ = 0) sites. 
Density matrix renormalization group diagonalization of finite systems (up to 120 sites) is used 
to calculate the charge-density correlation function and structure factor in the ground state. The 
modulation period can still be predicted through effective Fermi wavevectors, kp* , and densities, 
and we have found that it is much more sensitive to electron (or hole) doping, both because of 
the narrow range of densities needed to go from g* = to vr, but also due to sharp 2fcF*-4fcF* 
transitions; these features render CDW's more versatile for actual applications in heterostructures 
than in homogeneous systems. 

PACS numbers: PACS: 71.45.Lr, 73.21. Cd, 71.27.+a, 71.10.-w, 73.20.Mf, 72.15.Nj. 



The appearance of charge-density waves (CDW's) in 
conductors with a chain-like structure is a well known 
phenomenon) ''■■'?''j electrons condense into a ground state 
with a modulated charge distribution, as opposed to the 
uniform distribution of normal metals. Depending on 
their modulation period, these CDW's may slide along 
the chains to exhibit a wealth of interesting transport 
phenomena such as non-ohmic behavior, and AC current 
oscillations induced by DC electric fieldsiSi^ Nonethe- 
less, the fact that CDW transport was only controlled 
in bulk samples has hindered its integration to devices 
for some time. However, advances in material synthesis 
and manipulation over the last decade or so have given 
rise to new, micro- and nanostrucured, CDW systems. 
For instance, films of Rbo.soMoOa ("blue bronze") with 
thickness in the range 100-1000 nm have been grown by 
pulsed-laser deposition;^ mesoscopic wire structures have 
also been obtainedi^ Also, regions of negative absolute 
resistance have been observed at sufficiently low tem- 
peratures in single crystals of NbSea and o-TaSs with 
cross sections of 0.2 to 1 jjlw? .^■'' These approaches open 
up the possibility of fabricating CDW heterostructures, 
the properties of which are still unknown, and certainly 
worth investigating. 

While most of the above-mentioned features of CDW 
systems in confined geometries are dynamical, the knowl- 
edge of static properties is also of primary importance. 
From the theoretical point of view, the modulation pe- 
riod of the CDW is one of the key issues, even in the case 
of homogeneous systems. This has been addressed with 
the aid of various models, the simplest one being the one- 
dimensional Hubbard model, but the modulation period 
is still a matter of debate. On the one hand, the Luttinger 
liquid (LL) descriptior&SiiSiii predicts that charge corre- 
lations should be dominated by a wavevector q* — 2kp, 
where kp = TTp/2 [kp = 7r(2 — /o)/2] is the Fermi wave 
vector for an electronic density p < 1 [p > 1] of free 



electrons on a periodic lattice. By contrast, numerical 
analyses predict that the CDW instability is dominated 
by the Akp term, for strong enough on-site repulsion, 
jjji2ji2,i4ji^ These two scenarios can be reconciled if the 
amplitude of the 2kp contribution vanishes above some 
crossover coupling, U^{p). At any rate, considerable in- 
sight can be gained by examining the behavior of CDW's 
in a special class of heterostructures. Indeed, a simple 
model to study the interplay between electronic correla- 
tions and an underlying superlattice structure has been 
proposedjiS, its Hamiltonian reads 
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where, in standard notation, i runs over the sites of a 
one-dimensional lattice, c^^ {cia) creates (annihilates) a 
fermion at site i in the spin state a =\ or j, and riia = 
c^^Cia. The on-site Coulomb repulsion is taken to be 
site-dependent, in a repeated pattern of Ljj sites with 
Ui — U > Q (the repulsive layer), followed by Lq sites 
with Ui = {) (the free layer); one defines the 'aspect ratio' 
as a, — Lu/Lq. 

The magnetic properties of the above model turned 
out to be quite different from those of the corresponding 
homogeneous systemiiSjiiii^ Also, the density at which 
the system is a Mott insulator, pj, is shifted from half 
filling to a value which depends on the aspect ratioiiS 
These findings on a discrete lattice have been confirmed 
by an extension of the above model to a superlattice made 
up of Luttinger liquids (LL) ^QtSL 

As far as charge-density waves in these Hubbard super- 
lattices are concerned, a preliminary account of their be- 
havior has been presented in R,ef.l22t by virtue of the lim- 
itations imposed by Lanczos diagonalizations, that study 
was restricted to the electronic density p = 11/6 on lat- 
tices with Ns < 24 sites. Nonetheless, the results were 
quite surprising: it was found that the CDW wavevector 
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Figure 1: (Color online) The charge-density structure fac- 
tor, for two different system sizes, A'^s = 48 (circles) and 96 
(triangles), and for the same density: the cusp is located at 
g* = 7r/4, irrespective of Ns. The inset shows its Fourier 
transform, the correlation function, for the 48-site superlat- 
tice, with r = taken at a free site: the curve guides the eye 
through the oscillatory part, whose period is marked by the 
double arrow. 



can be cast in a simple form, q* = 4A:f*, where the ef- 
fective Fermi wavevector, kp* , depends on the density p, 
on Lf/, and on Lq (see below); also, this CDW wavevec- 
tor oscillates with the number of spacer (i.e., free) sites, 
with a period which is half of that for the spm-density- 
wave for the same case. In view of this, many interesting 
questions arise, such as how does q* depend on the elec- 
tronic density and on the Coulomb repulsion U, and if 
there are sharp 2/ci?*-4/ci?* transitions. In order to an- 
swer these questions one needs to examine both a wide 
range of band fillings and large lattices, so here we report 
on a Density Matrix Renormalization Group (DMRG) 
approach22i2iiSSi2& to this problem. 

We consider the Hamiltonian ^ on lattices with Ns 
sites, Ne electrons, and open boundary conditions. The 
appropriate Finite Size Scaling (FSS) parameter, how- 
ever, is the number of periodic cells, Nc = Ns/Ni,, for 
a lattice basis with Ni, = Lu + Lq sites. The ground- 
state, iV'o), and energy are obtained through the DMRG 
method. Lattice sizes ranged from 48 to 120 sites and 
truncation errors in the DMRG procedure were kept 
around 10~^ or smaller. We have performed a systematic 
study of the charge-density structure factor. 



ciq) = ^T.^''^''' 



\QrQj) 



(2) 



*j 



where {QiQj) = {'(po\ninj\ipo) withrii = ni-^+riii; the sig- 
nature of a CDW instabihty is a cusp in C{q) at q = q* . 
(We chose to examine the charge-density structure fac- 
tor without the uncorrelated part removed; in so do- 




Figure 2: (Color online) The charge-density structure factor 
at different electron densities, for a superlattice with 24 cells 
(48 sites): in each case, filled and empty symbols respectively 
mark 2k f* and Akp* (see text), and each plot has been shifted 
vertically by the same amount of 0.1, for clarity. 



ing, C(0) = N^/Nc, instead of satisfying the sum rule 
C(0) = 0, but we benefit from the fact that the cusps are 
much sharper, and thus more readily located.) We have 
considered different values of the Coulomb repulsion U, 
different occupations p — N^/Ng and different configura- 
tions {Ui}. Not all configurations {Ui} fit into all sizes 
and occupations but, since DMRG allows us to study 
a wide range of lattice sizes, we were able to establish 
overall trends. 

A typical plot of C{q) is shown in Fig. 2] for lattices 
with 48 and 96 sites, Ljj = 1, Lq = 2, and for density 
5/4: a cusp is observed at q* = tt/A. It is also evident 
from the figure that neither the cusp position nor the cusp 
height depend on the system size; thus, our analyses are 
not plagued by subtle finite-size effects. The inset of Fig. 
^ displays the charge-density correlation function in the 
case of Ns =48; only the correlation between charges on 
free sites are plotted, in order to separate intercell from 
intracell oscillations. We see that once boundary effects 
die out (typically beyond three unit cells) an oscillatory 
regime sets in with period A* = 27r/q* = 8 unit cells. 

Before comparing the behavior at different electron 
densities, we recall some of its special values which play 
a crucial role in establishing different regimesiiSiSS First, 
Pq = l/{Lij + Lo) signals the presence of one electron per 
cell. For U large enough, added electrons will accumu- 
late on the free sites; this saturates at p-^ = 2/(1 -I- i'). 
As the density increases further, one eventually reaches 
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corresponding to the addition of one electron to every 
repulsive site. Also, the cell electronic density is 
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Figure 3: (Color online) The cusp position as a function of the 
electrinc density, for three representative superlattice config- 
urations, all with U — 4: (a) for Lu = Lo = 1, (b) for 
1/(7 = 1, Lo = 2, and for Lu = 2, Lq = 1. The points have 
been determined through analyses of C{q), while the dashed 
and full straight lines respectively represent 2kp and ikp [see 
Eq. ®]. 
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Figure 4: (Color online) The charge-density structure factor 
at a fixed electron density (p = 7/4), and for several values 
of U, for a superlattice with 32 cells (96 sites): in each case, 
filled and empty symbols respectively mark 2fc_F* and Akp* ■ 



Figure 121 shows C{q) for the Ljj = Lq = 1 superlattice 
for a wide range of densities, all with U = 4. The cusps 
move as the density varies, but two regimes can be clearly 
distinguished: for p < p^^ = 1, the cusps occur at q* = 
2k F*, whereas for p > p^ = 1, they are at q* = 4k p*, 
where the effective Fermi wave vector, kp, is defined as 
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vrpccii if Pccii < 1 

7r(2 - pccii) if Pcoii > 1- 
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The above analysis was repeated for different super- 
lattice configurations, and we followed the evolution of 
the cusp positions with the electron density. Some of the 
results (still for U = 4) are summarized in the q* vs. p 
plots shown in Fig. |3| The signature of the 2kp*-4kF* 
transition is now a doubling of the slope. For this value 
of U and for Ljj < Lq, the transition occurs a,t p = p-\i, 
which can be understood by noticing that electrons oc- 
cupy the free sites for densities below p^; above this 
density, additional electrons occupy the repulsive sites, 
thus forcing an overall charge rearrangement in order to 
minimize the on-site repulsion. For Ljj > Lq, on the 
other hand, p|| < 1 so that electrons can occupy repul- 
sive sites at no cost in energy until one reaches p — pi; 
in this case, it is only for densities larger than pi that a 
major charge distribution becomes imperative. Another 
important feature of a superlattice structure is that all 
modulations become accessible within a range of dopings 
smaller than in homogeneous systems, for which q* = irp 
or 2ttp, corresponding to 2fcj? or 4kF, respectively. 

While the scenario of Fig.Ois generic for all U > 4, the 
same does not hold for [/ < 3. Figure 0] shows C{q) for 
several values of f/, at a fixed density p > p|| : we see that 
the cusps move from 2kp* to 4kF* as U is increased. As 



the density is further increased, the 2kF*- 4kF* transition 
occurs at smaller values of C/; indeed, at weak coupling 
the occupation of repulsive sites can be compensated by 
delocafizing the electrons. 

The picture that emerges is then summarized in Fig. 
O with the exception for the case Ljj = Lq = I, for all 
other superlattices the transition continuously migrates 
to larger densities, along a phase boundary line, Uc{p)- It 
should be noted that the transition boundary for p < pi, 
in the case of L[/ > Lq, is not as steep as that for Lj/ < Lq 
and p < p||; for intance, for the specific case of Lij = 2 
and Lq = 1, we estimate Uc ~ 20, 8, and 3, for p = 7/6, 
5/4, and 17/12, respectively. 

For consistency, we have also examined the Luttinger 
liquid (LL)iiiSi critical exponent, Kp, which is relatedSS 
to C{q) through 
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From a plot similar to that of Fig.[21(i.e., one for the same 
parameters, but with the uncorrelated part removed) , we 
see that the slope at g = first increases with p, reaches 
a maximum value and then decreases to zero at p/; the 
actual values of Kp are in excellent agreement with those 
obtained in the LL approach to this problem. ^"-^^ By the 
same token, we have found that Kp steadily decreases 
with U at fixed density, similarly to the homogeneous 
LL. 

In summary, the present study has established that 
superlattices made up of normal metals and of charge- 
density-wave materials lead to a wider choice of outcomes 
for the CDW period than in homogeneous systems: the 
modulation period is much more sensitive to electron (or 
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Figure 5: (Color online) Schematic phase diagram for the 
2k F*- 'ikp* transition in the space U vs. p: The top two panels 
correspond to superlattices with the same aspect ratio, £ = 1, 
while the bottom two are for £ < 1 (left) and £ > 1 (right). 



dence with the homogeneous systems. Indeed, between 
intermediate and strong couplings, these transitions are 
density driven, and occur when the electronic density is 
increased beyond one electron per free site, if there are 
more free sites than repulsive ones; if there are more re- 
pulsive sites, the transitions take place near the Mott- 
Hubbard insulating density. For weak couplings, and 
if the number of repulsive and free sites are different, 
the superlattice structure plays a secondary role in de- 
termining the modulation: the critical on-site repulsion 
displays just a weak dependence with the density; in this 
regime, the 2kp*-Akp* transition may be driven by ap- 
plying pressure, which increases U. Therefore, from the 
standpoint of applications, CDW's seem to be more ver- 
satile in heterostructures than in homogeneous systems, 
and experimental work aiming at assembling CDW ma- 
terials along these lines should unveil a variety of new 
effects and properties. 
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